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ABSTRACT 

We find that a steady state bar induces transient features at low velocities in the solar neighborhood 
velocity distribution due to the initial response of the disc, following the formation of the bar. We 
associate these velocity streams with two quasi-periodic orbital families, librating around the stable 
a;i(l) and xi{2) orbits near the bar's outer Lindblad resonance (OLR). In a reference frame moving 
with the bar, these otherwise stationary orbits precess on a timescale dependent on the strength of the 
bar, consistent with predictions from a simple Hamiltonian model for the resonance. This behavior 
allows the two orbital families to reach the solar neighborhood and manifest themselves as clumps in 
the u-v plane moving away from (xi(2)), and toward (a;i(l)) the Galactic center. Depending on the 
bar parameters and time since its formation, this model is consistent with the Pleiades and Coma 
Berenices, or Pleiades and Sirius moving groups seen in the Hipparcos stellar velocity distribution, if 



the Milky Way bar angle is 30 



° < 



< 45° and its pattern speed is rifc/rio = 1-82 ± 0.07, where fio is 



the angular velocity of the local standard of rest (LSR). Since the process is recurrent, we can achieve 
a good match about every six LSR rotations. However, to be consistent with the fraction of stars in 
the Pleiades, we estimate that the Milky Way bar formed ~ 2 Gyr ago. This model argues against a 
common dynamical origin for the Hyades and Pleiades moving groups. 
Subject headings: stellar dynamics. Galactic bar, solar neighborhood 



1. INTRODUCTION 

The Milky Way bar is hard to observe directly due to 
our position in the Galactic disc. Thus, its parameters, 
such as orientation of its major axis with respect to the 
Sun-Galactocentric line and angular velocity, have been 
inferr ed indirectly from observations of the in ner Galaxy 
(e.g., IBlitz and SpCTgellTMt IWeinberdfTOOl) . 

However, numerical modeling suggests that the bar 
can also affect the local stellar velocity distribution. 
It has been found to account for the Hercules stream 
(jCJ,uillcn 2003; Fux 2001; Dehnen 2001^, the vertex devi- 
ation (jMii hlbaucr and Dehnen 2 0031) . errors in the Oort 
constants (.Minchev et al.. ,2007i) . and (in combination 
with spiral structure ) the observed flatness in t he age- 
metallicity relation ([Minchev and FamaevI [20091 ) . The 
bar angular velocity, or pattern speed, Ob, has been 
well established to be such that the solar circle lies just 
outside the ba r's 2:1 outer Lindblad resonance (OLR). 
iDehnenI ([20000 found that in order to reproduce the Her- 
cules stream, the bar pattern speed should be rib/Oo = 
1.85 ± 0.15, where VLq is the angular velocity of the 
local standard of rest (LSR), and the bar orientation 
with respect to the Sun-Galactocentric line is in the 
range 10° < < 70°. B y accounting for tren ds seen 
on the Oort constant C, iMinchev et al.l ([2007D found 
f]b/flo = 1.87 ± 0.04 and 20° < 4>o < 45°. Such bar 
orientation is consistent with estimates derived from IR 
photometry (15° < < 45°) and OGLE-II microlens- 
ing observations of red clump giants in th e Galactic bulge 
(24° < 00 < 27°. lRattenburv et al.|[2007[) . 

Recent analysis of the local velocity field show that 
the origin of the most prominent low-velocity moving 
groups in the solar neighborhood favor a dynamical ori- 
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gin ([Famaev et al.l[2008h . While the bar has been shown 
to affect the u-v plane at higher velocities, spiral den- 
sity waves have been found to create resonant structure 
in the lower velocity regions, such as the splitting of 
the Pleiades /Hyades and Co ma Berenices moving groups 
([Quillen and MinchevI 120051 ) . Is it re ally t rue th at the 
bar is not effective at low velocities? iFuxl ([2001[ ) inves- 
tigated the effect of the bar on the local veloci ty distri- 
bution, similarly to the work by IDehnenI ([2OOO0 . but for 
test particles integrated forward in time. The author re- 
ports time varying structure in the u-v plane following 
the growth of the bar. However, he did not investigate 
this further, being more interested in the effect of the bar 
once the disc has fully responded to the perturbation. To 
increase the particle statistics and reduce ph ase mixing , 
the u-v distributions presented in the work bv lFuxl ([20011 ) 
were averaged over 10 bar rotations and then smoothed. 

However, it is possible that the effect of the initial 
response of the disc to the formation of a central bar 
can still be seen today in the solar neighborhood veloc- 
ity distribution, provided the bar has formed or evolved 
recently. There are both observational and theoreti- 
cal arguments suggesting such a scenario. By using 
data fr om the Two Micron All-Sky Survey (2MASS), 
iCole a nd Weinberg (2002) estimated that the Milky Way 
bar is likely to have formed more recently than 3 Gyr 
ago and suggested that this event could have been tri; 
gered by a now-merged satellite. IMinchev et al.l ([200' 
showed that phase wrapping in the thick disc caused by 
initial conditions which might have been left following a 
merger, can explain four observed high velocity streams. 
To match the locations of these streams, a strong per- 
turbation of the MW disc is required about 2 Gyr ago. 
Such an event could have triggered the formation of the 
MW bar. In addition, cosmological simulations show 
that massive minor mergers are likely to have happened 
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Fig. 1. — Time development of the u — v plane for a simulation of a stellar disc with Milky Way kinematics perturbed by a steady state 
bar. The pattern speed is fixed at Qi,/Qo = 1.85, the maximum sample depth is dmax = 250 pc, and the bar orientation is 0o = 25° 
(geometry as shown in figure [2]l . Shaded contours show the particle number density. Each panel is a snapshot of a particular time of the 
same simulation, in increments of one rotation at t-q. Note that even though the bar perturbation is stationary, structure varies periodically 
with time. 



during the lifetime of M i lky W ay size system. For ex- 
ample, |Kaia^iidii^^F^L| ()2008D estimated that as many 
as five objects more massiye than 20% Mdisc could have 
been accreted since 2; ~ 1. All this evidence from obser- 
vations and simulations calls for a detailed investigation 
of the effect of a young, or a recently evolved bar, on the 
solar neighborhood velocity distribution. 

In this paper we examine the effect of a central bar on 
the local velocity distribution by simulating the time evo- 
lution of a barred stellar disc. Unlike in previous works, 
here we do not time average over particles positions, in 
order to investigate in detail the response of the disc. 
We explore the possibility that stellar streams in the u- 
V plane could be induced by a recent bar formation or 
evolution. 

2. SIMULATION SET UP 

We perform test-particle simulations of a stellar galac- 
tic disc with parameters consistent with Milky Way 
kinematics. The density distribution is exponential, 
S(r) ^ e~''/''f, with a scale length = 0.38rn , con- 
sistent with the 2MASS photometry (Oihal 120011 ). We 
give particles an initial radial velocity dispersion in the 
form of a Gaussian distribution. The Milky Way disc 
is known to have a radial velocity dispersion which de- 
creases roughly exponentially outwards: cr^ ~ . 
In accordance with this we implement an exponential 
decrease in the standard deviation of the radial veloc- 
ity dispersion of stars with radius, with a scale length 
r„2 — 0.45ro (|Lewis and FreemanI Il989[ ). Since we are 
interested in the low-velocity regions of the u-v plane, 
we set (T„ = 10 km/s at tq. 

To explore the time development of the system, we do 
not time-average over position and velocity vectors, as it 
is frequently done in. test- particle simulations (e.g.. lFru3 
I2OOII : iMinchev et al.|[2007[ ) where no dynamical develop- 
ment of the system is expected. For each simulation we 



integrate N=5 x 10^ particles for 30 rotation at the solar 
radius, rg. 

The background axisymmetric potential due to the disc 
and halo has the form $o(^) = ^olog(r), corresponding 
to a flat rotation curve. We model the nonaxisymmetric 
potential perturbation due to the Galactic bar as a pure 
quadrupole 



= Ab(eb)cos[2(0-f7b<)] X 




(1) 

Here Ab(eb) is the bar's gravitational potential ampli- 
tude, identica l to the same name parameter used by 
iDehnenI ()2000[) : the strength is specified hy t\, — —a from 
the same paper. The bar length is rb = O.Srcr with Tcr 
the bar corotation radius. The pattern speed, i^b is kept 
constant. The bar amplitude e is initially zero, grows lin- 
early with time at < i < ti and transitions smoothly 
to a constant value after t ^ ti = A bar rotations. This 
insures a smooth transition from the axisymmetric to the 
perturbed state. 

In our units the solar neighborhood (SN) radius is rg = 
1; the circular speed is = 1 everywhere since the rota- 
tional curve is flat. To convert to real units we use Local 
Standard of Rest (LSR) tangential velocity of 240 km/s, 
and Galactocentric distance of 8 kpc. The 2:1 OLR with 
the bar is achieved when fib/^^o = 1 + k/2 « 1.7, where 
K, is the epicyclic frequency. For a flat rotation curve 

3. RESULTS 

For an axisymmetric disc the potential is independent 
of time, thus energy and angular momentum are con- 
served. In the case of a single periodic perturbation, 
which is stationary with time, there is still a conserved 
quantity in the reference frame rotating with the pat- 
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Fig. 2. — Time development of nine, initially circular rings near the bar 2:1 OLR, in the bar reference frame. All simulation parameters 
are the same as in figure [T] The 2:1 ILR resonance and corotation are shown as dotted circles; the 2:1 OLR location is indicated by the 
dashed circle. The solar neighborhood is shown by a black dot, just outside the 2:1 OLR at an angle of 25° lagging the bar. Two families 
of orbits associated with the xi(l) (red) and xi(2) (blue) orbits process at different rates, giving rise to the streams in figure[l] 



tern. This is the Jacobi's integral, J = E — Lfli,, where 
E is the energy of the particle, L is its angular momen- 
tum, and rib is the pattern angular velocity. In such a 
system no ti me variations are expe c ted in phase space. 
For example, iMinchev and QuillenI (|2006( ) showed that 
no increase in velocity dispersion with time (disc heat- 
ing) occurs once the pattern has been fully grown. Con- 
sequently, in previous such studies a time-averaging pro- 
cedure of the particle positions has been employed in 
order to ach ieve better statistics and minimize co r nputa - 
tional time (|Fuxll2001l : lMinchev and Quillenll2007ll2008l ). 
Unlike in those works, here we look at the time evolution 
of the system by integrating a large number of particles 
and follow their evolution in phase space. 

3.1. Time evolution of the u — v plane 

Figure [T] shows the time development of the u — v plane 
for a simulation of a stellar disc with Milky Way kine- 
matics perturbed by a steady state bar. The pattern 
speed is flt/flo — 1.85, the maximum sample depth is 



dmax = 250 pc, and the bar orientation is 0o = 25° (ge- 
ometry as shown in figure [J). Shaded contours show the 
particle number density N . Each panel is a snapshot of 
a particular time of the same simulation, in increments 
of rotation at ro, up to < — 23. 

The bar is fully grown in four bar rotations, which 
corresponds to t « 2.2 in the figure. Note that even 
though the perturbation is stationary, structure varies 
with time. At i « 6 the distribution becomes bimodal, 
forming a strong feature at {u,v) ~ (—10, —20). During 
the next five rotations this clump moves toward more 
negative radial velocities, becomes weaker and eventu- 
ally disappears at about (u,t) — (40, 11). This process is 
repeated approximately every 6.5 rotations, starting at 
times t sa 6, 12 and 19. In addition, there is a second fea- 
ture forming at t « 7, which evolves with time to more 
positive radial velocity m, and approximately constant 
V ~ —10 km/s. This clump recurs on a slightly shorter 
timescale (At w 5.5) than the one described above, as ap- 
parent from the two slopes formed by aligning the clumps 
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Fig. 3. — Top: Change in average radius with time for the ini- 
tially circular rings shown in figure [2] The solar and OLR radii 
are indicated by the dashed lines at r = ro and r 0.92ro, re- 
spectively. Bottom: Same as the top panel but for two annuli of 
radial width Ar = 0.02ro just inside (blue, solid) and just outside 
(red, dashed) the 2:1 OLR. Red and blue orbits correspond to the 
streams at m > and n < in figure [l] respectively. The ampli- 
tude of the oscillations decreases with time as the system relaxes, 
similarly to the strength of the velocity streams in figure [T] 

at w < and m > 0, for example in the fifth column in 
figure [TJ One should also note that the strength of the 
features diminishes at every period. 

Note that the structure in velocity space seen in figure 
[1] is not due to our ICs as we initially integrate particles 
in the axisymmetric potential for 3 Gyr to insure the disc 
is relaxed before the bar is grown. The reason for this 
transient nature of velocity streams in the u-v plane is 
an effect caused by the initial response of the disc to the 
recent formation of the bar. To show this we next look 
at the time evolution of orbits near the 2:1 OLR. 

3.2. Orbits near the 2:1 OLR of the bar 

Linear theory predicts that the orientation of orbits 
shifts across the 2:1 OLR of the bar, from perpendic- 
ular to the bar major axis insid e the OLR, to paral- 
lel w ith it outside the resonance (|Binnev and Tremaind 
I2OOI . These are referr ed to as the xi ( 2) and xi ( 1 ) orbits, 
respe ctively (see, e.g., lQuillen|[2003HFia^ 120011: iDehnenI 
I2OOOD . Near th e peak of the resoiiance, both typ es of 
orbits can exist fContopoulosl llQTSl: lWeinberdri994[ ). As 
suggested by Kalnais (1991), for a SN location near the 
OLR and a bar angle in the range 0° < 0o < 90°, the 
closed orbits from either side of the OLR would pro- 
duce two streams; one moving inward {u > 0) and the 
other outward {u < 0), which could be associated with 
he Hyades and Sirius stellar streams seen in the Hippar- 
cos stellar velocity distribution. Is it possible that the 
two streams we observed in figure [T] above are somehow 
linked to these orbital families? To try to answer this 
question, let us look at how orbits evolve with time. 

In figure [5] we show the time development of nine ini- 
tially circular rings near the bar 2:1 OLR, in the bar ref- 
erence frame. For each ring particles start at the same 



initial galactic radius and with a random distributed in 
azimuth. All simulation parameters are the same as in 
figure [1] Rotation is clockwise, thus in the reference 
frame of the bar particles move in the counterclockwise 
direction. The bar is represented by the ellipse in the cen- 
ter of each panel. The location of the 2:1 ILR and the CR 
are shown as dotted circles and the 2:1 OLR location is 
indicated by the dashed circle. The solar neighborhood 
is shown by a black dot, just outside the 2:1 OLR, for 
a bar angle of 25° (as in figure [T]). The initial radii are 
such that four rings lie just inside the OLR (blue end) 
and four are just outside the resonance (red end). One 
ring is situated right on the resonance (green). 

As the system develops with time, the inner and outer 
groups of rings tend to stay individually aligned and 
precess at different rates. Since these also exhibit ra- 
dial oscillations, in order to see the timescale better we 
plotted the average radius of each ring as a function of 
time in figure [31 top panel. What we would actually 
see near the Sun is the average contribution from these 
orbits. Expecting that the ones just inside the OLR be- 
have differently than those just outside the resonance, 
we want to know what is the frequency of oscillation 
for each group. By considering annuli of finite width 
Ar = 0.02ro (bottom panel of figure [31), we can now see 
that the timescales of precession inside and outside the 
OLR are about 5.5 and 6.5, respectively. These are con- 
sistent with the recurrence times of the velocity streams 
in figure [H Thus, we can associate the stream at w < 
with a quasi-periodic family of orbits librating around 
the xi(2) closed orbit. Similarly, the stream at u > 
originates from an orbital family librating around the 
xi(l) close orbit. 

The amplitudes of oscillation inside and outside the 
OLR seen in the bottom panel of figure [3l decay with 
time; in other words, the disc relaxes. This is consis- 
tent with the decrease in the fraction of particles in the 
streams in figure [H as each cycle is repeated (see, e.g., 
t = 7,13 and 20). However, stellar streams in the u-v 
plane remain quite strong until the end of the simulation 
at i = 23 LSR rotations. This shows that the disc re- 
laxation time is much longer than the age of the Galaxy. 
We also expect to see the effect of phase wrapping. For 
example, at times t > 17 in figure [1] two steams are ap- 
parent at w < where the same orbital family is seen at 
the end and beginning of its precession around the xi(2) 
orbit. However, in general it appears th at phase wrap- 
ping e ffects are not important, unlike in iMinchev et al\ 

3.3. Variation with bar pattern speed and orientation 

Figure [4l shows simulated SN velocity distributions 
similar to those shown in figure [1] Here however, while 
the contours still show the particle number density, the 
color levels show the change in angular momentum AL = 
L — Lq, where L and are the current and initial angu- 
lar momenta. The AL color levels are normalized for all 
panels, allowing one to see the relative change in angular 
momentum for the full range of parameters shown. The 
color bar values can be converted to (pc km/s) by mul- 
tiplying by IOOvq. From top to bottom 0o changes from 
to 160° in increments of 20°. Different columns show 
different pattern speeds ranging from fib/r^o = 1.7—1.95 
(form left to right). The solar radius is at the bar's OLR 
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Fig. 4. — Variation in the u — v plane with bar pattern speed and orientation. Contours show particle number density, while the color 
levels represent the change in angular momentum AL for a given location on the u-v plane. Color bar values can be converted to (km/s 
pc) by multiplying by lOOtio. Different rows show different bar angles from ij>0 = 0° to = 160° in increments of 20°. Different columns 
show changes in bar angular velocity for solar circle closer to the 2:1 OLR (fi^ = 1.75f!o) to farther away from it (fij = 2.0f!o). 



in the leftmost column. Particles are confined to a region 
of 250 pc from the Sun. 

It is important to realize that for slower bars (SN closer 
to the resonance) structure develops later than for faster 
ones (see figure [5]) . Consequently, snapshots of the u-v 
plane at the same time since bar formation and different 
angular velocities will show different stages of the preces- 
sion cycle. To approximately account for this, in figure U 
we show snapshots at earlier times as the pattern speed 
increases. For example, when i7f, = 1.7Qq we show a plot 
at t = 10, while for fit, = IMflo t = 5. 

Note the strong variation with both pattern speed and 
orientation. We first look at the variation with angle for 
the same simulation as in figure [1] shown in the forth 
column of figure The clump at u < which we asso- 
ciated with the orbital family just inside the OLR (see 



section [3. 2p . remains at a constant tangential velocity v 
as (f)o changes, but shifts in the positive u direction. It 
is intriguing that at angles 4>o > 60° this stream moves 
in the positive radial velocities region. At the same time 
the stream associated with the a:i(l) closed orbit shifts to 
negative u. A similar behavior is observed for all other 
pattern speeds as well. This is easy to understand in 
view of the precession (in the reference frame of the bar) 
and pulsation of the orbital families just inside and just 
outside the OLR (see figure [2]). By changing the SN an- 
gle, the two families of orbits reach the solar vicinity at 
different orientations, which results in sampling different 
subsets of the orbits in the u-v plane. 

Examining a particular row in figure SI we see that as 
Qb increases the streams shift to more negative angular 
velocities. This is expected if the stars in the streams 
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formation ti in units of bar rotations. The y-axis shows the time 
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come from near the OLR radius, since the distance from 
the solar circle to the OLR circle increases with an in- 
crease in We can see this also in the change in an- 
gular momentum shown by the color levels in figure 01 
The largest AL occurs for $7^ = l-9^o, meaning that the 
guiding radius for these stars is smallest. Note that since 
the expected bar angle is in the range 10 < 0f, < 60, pos- 
sible matches to the observed velocity field will always 
have It < for the xi(2) family and u > for the a;i(l) 
one. 

As apparent form examining the fifth row in figure S] 
we see that beyond fib = LSSfio phase wrapping is more 
efficient near the resonance and becomes unimportant 



away from it. 

3.4. Variation with bar growth time 

In this paper in all figures we grow the bar in four bar 
rotation. In order to see how our results would change 
with a change in the bar growth time ti, we ran simu- 
lations with different values in the range < < 12 
in units of bar rotation. We plot the results in figure [5j 
The y-axis shows trmax , the time for which the functions 
in figure |3] peak for the first time. The x-axis shows the 
bar growth time used for that particular simulation. Red 
and blue symbols correspond to the orbital families just 
outside and just inside the bar's OLR. The dependence is 
almost linear with trmax ~ l/4ti-|-3. The relative phases 
of the functions in figure [3] remain the same for all values 
of ti, as do all our other results. We conclude that while 
this delay in structure formation as ti increases would 
introduce uncertainty in our prediction for the bar for- 
mation/evolution time, the general effect of bar growth 
time in unimportant. 

3.5. Variation with bar strength 

Our default value for the amplitude of the bar per- 
turbation is e = —0.012. We would like to know how 
changing this will affect our results. We first derive an 
analytical expression for the orbital libration timescale 
At, as a function of bar strength e. Next we compare 
t hat to our results from simulations. 

iContopou los (1975, 1988) showed that the dynamics 
of stars confined to the Galactic plane in the case of an 
axisymmetric potential can be described by a Hamilto- 
nian written in a thir d order post-ep icyclic approxima- 
tion. Taking it further. lOuillenl ()2003t ) considered the dy- 
namics of stars that are affected by perturbations from 
both spiral structure and the Milky Way bar. She con- 
structed a simple one dimensional Hamiltonian model 
for the strongest resonances in the epicyclic action angle 
variables. We can use her equation 22 to describe the 
effect of the bar only, by neglecting the term related to 
spiral structure (last term on the right hand side). In a 
slightly different notation we can write 

H = ap^ + Sp + Pp^I'^ cos 0, (2) 

where p is the action variable related to eccentricity and 
(j) is the r esonant angle . The coefficient a = — 5.7r(^^ 
(eq. 18 in 'Qui lle"nll2003| ) is constant at the solar radius 
and (S = K + 2(ri — f^h) sets the distance from resonance. 

For a flat rotation curve (3 « QMQ^j2fil{Af /ro){rblroY 
and is related to the bar amplitude e through Af = 
(ef;^/3)(ro/rf,)3 (eq. 7 in IDehnenI f2000f) . where r^ ~ 
0.43ro is the radius at which the bar ends (for Vlb — 
1.85rio). By dimensional analysis from a and j3 we can 
derive a timescale for the resonance: 

Ai« (27r)-ia-i/3/3-2/3^ (3) 

where we have divided by 27r to put the expression in 
units of LSR rotations. This is the libration time around 
the fixed points associated with the a;i(l) and j;i(2) or- 
bits. Note that the functional dependence of libration 
time on bar strength is 

At ~ t~^/\ (4) 
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Fig. 7. — Similar to figure [T] but showing variations with angle (left to right) and time (top to bottom) for the same pattern speed 
f2(, = l.SSHo. For this pattern speed, possible matches to the Pleiades and Coma Berenices streams are achieved for t = 7, 8 and 
35° < 00 < 45°. 
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Fig. 8. — Top left: Hipparcos stellar velocity distribution 
with the Sun's motion subtracted (values from lDehnen and BinnevI 
I199SI 1. Top right: Simulation with initial velocity dispersion 
CTt- = 40 km/s and parameters as indicated in the figure. Good 
match to Sirius and Pleiades groups. Bottom right: Small changes 
in the bar's orientation and time since formation, provides good 
match to the Coma Berenices and Pleiades groups. Bottom left: 
As in figure [2] with simulation parameters as the corresponding 
u-v plot on the right. Arrows show the orbital families giving rise 
to the streams in the velocity distribution. 

To compare to this analytical result we performed sim- 
ulations with different perturbation strengths, in addi- 
tion to the default value of e = —0.012. In figure [B] the 
open blue squares and the red star symbols show the 
precession time for the streams just inside and just out- 
side the OLR for a bar pattern speed as in figures [T]- 
[3] and different bar perturbation amplitude in the range 
-0.006 ^ e ^ -0.018. The sohd line plots equation 
[3] scaled by a factor of 0.38 (estimates by iQuillen ,2003i 



are accurate to within an order of magnitude). We see 
excellent agreement in the functional behavior between 
our numerical and analytical results. Simulations with 
stronger bars result in shorter libration time. Features 
in the u-v plane remained of similar strength to those 
seen in figure [TJ 

4. CONSTRAINING THE BAR 

As we showed in section [3?2| for a pattern speed consis- 
tent with expectation for the MW bar, we can associate 
two streams in velocity space with the quasi-periodic 
families of orbits librating around the xi(2) and a;i(l) 
closed orbits. In section [3731 we explored the variation of 
structure in the u-v plane with angular velocity and bar 
angle for a particular time. Now we would like to see 
how we can use this information to put constraints on 
bar parameters, such as pattern speed, orientation and 
formation time. 

For a given time, the tangential velocity v of resonant 
features in the u-v plane is set by the bar pattern speed 
as evident from figure HI On the other hand, the ra- 
dial velocity u, is set by the bar's orientation. Thus, 
assuming some features at low velocities in the Hippar- 
cos velocity distribution are of resonant origin, we can 
match stream positions in the u-v plane and estimate fib 
and (jjQ. In addition, structure varies with time due to 
the libration of the quasi-periodic orbits around the fixed 
points; this allows us to constrain the bar formation time. 
We need to compare our synthetic velocity distributions 
to the Hipparcos u-v plane corrected for the mot ion of 
the Su n. As recently discussed by McMillan and BinnevI 
(|2010). the velocity of the Sun with respect to the LSR, 
velocity of the LSR and the Sun's Galactocentric radius 
may be quite different than the currently accepted val- 
ues. However, the LSR angular rotation rate = vo/rg 
is found to be well constrained in the range 29.9-31.6 
km/s/kpc. r^o sets the relative position of streams in the 
u-v plane. Therefore, as we consider possible matches 
to the observed velocity distribution we need to be most 
concerned with the relative orientation of streams. In 
this paper we use the values vq — 240 km/s and tq = 8 
kpc resulting in fig = 30 km/s/kpc. 

In figure [7] we show how structure in the u-v plane 
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changes with bar angle and time, for the simulation in 
figure[T](Ob = 1.85f7o)- While the pattern speed is kept 
fixed everywhere, different columns show variation with 
angle in the range 25° < 4>o < 50° in increments of 5°. 
Rows show changes with time in units of rotations at tq 
for t = 6,7, and 8. The clump at negative u can be asso- 
ciated with the Pleiades moving group, while the one at 
u > at the beginning of its formation is consistent with 
the Coma Berenices stream. For this pattern speed, good 
matches are achieved for t = 7,8 and 35° ^ (j)Q S 45°. 
Since structure in the u-v plane is recurrent about every 
six LSR rotations (see figure [ij, one would expect later 
times to be equally consistent. However, we require a suf- 
ficiently strong signal in order to account for the fraction 
of stars in the Pleiades, thus we reject later times. 

In figure |8] we compare the observed velocity dis- 
tribution to our simulation s. We used the catalog 
from Holmb erg et al.l (|2009f ) to plot the u-v distribu- 
tion of Hipparcos stars in the top left panel of figure 
11 To subtract the S olar m otion we used the values by 
Dehnen and BinnevI ()1998[ ). In the top right panel we 
show a simulation with parameters as indicated in the 
figure. Note that here we start with a radial velocity 
dispersion cr„ = 40 km/s in order to populate the Her- 
cules stream at v «= 40 km/s. The two clumps in the 
top right panel are consistent with the Pleiades and Sir- 
ius groups. For a 10° change in bar orientation and half 
a rotation at ro we can associate these groups of stars 
with the Pleiades and Coma Berenices streams (bottom 
right). Finally, the bottom left panel in figure [S] shows 
the orbital families giving rise to the streams in the u-v 
plane. 

Similar inspection of the u-v plane for pattern speeds 
in the range 1.7 ^ ilb/ilo ^ 2.0 reveal that, depending 
on the bar parameters and time since bar formation, we 
can obtain a match to either the Coma Berenices or Sir- 
ius moving groups in addition to the Pleiades stream. 
We conclude to be consistent with structure in the ob- 
served u-v velocity distribution our model requires the 
bar angle to be in the range 30° ^ (j>o ^ 45° and its 
pattern speed to be flb/flQ = 1.82 ± 0.07. This is in very 
good agreement with previous estimates of these bar pa- 
rameters. For example, by a quantitative comparison 
of the observ ed with the simulated position of the Her- 
cules stream, iDehnenI ([2000) deduced the MW bar pat- 
tern speed to be r2(,/r2o — 1-85 ± 0.15 and the bar angle 
10° ^ 00 ^ 70°. By relating the dynamical effect of the 
bar to the derivatives of the velocities via the Oort con- 
stant C, an additional cons t raint on bar parameters was 
provided bv lMinchev et al.l ()2007D . where they estimated 
flh/no = 1.87 ± 0.04 and 20° < 0o ^ 4 5°. Our result is 
also in agreement with the estimate bv iDebattista et aD 
()2002t ). based on OH/IR star kinematics. 

To account for the strength of the Pleiades, we find 
that the bar formation time has to be between 7 and 9 
rotations at tq, which corresponds to 1.75-2.25 Gyr. 

4.1. Relating to other work 

Previous works have also reported features (or not re- 
ported but are identifiable in their figures) in their barred 
disc simulations, which we can now explain as the initial 
response of the disc to bar for mation. 

One example is the work by 'Bagl ev et al.l ()2009D who 
investigated the effect of a central bar on disc morphol- 



ogy. In their figure 3 they present number density plots 
of the time evolution of a barred galactic disc, similar 
to our figure [2l but showing the extent of the full disc. 
In addition to the asymmetries in the R2 ring (associ- 
ated with the xi(2) orbit) reported by the authors, it is 
easy to identify also a second ring at smaller radius and 
different orientation relative to the outer one (a:i(l) or- 
bit). Even at the last three time outputs (figure ends at 
t = 12.5 in our units) where the outer ring appears sta- 
ble, close inspection of the region between the two rings 
reveals misaligned features similar to those seen in our 
figure [2] 

A second example is the work by iFuxl (|2001f) , where 
the author investigated the effect of the bar on the u-v 
plane. The second row of his figure 12 presents simulated 
local velocity distributions for a ratio of solar to OLR ra- 
dius ro/roLR. = 1-1 (or equivalently, fli, = 1.87f2o) and 
different bar angles. Due to the time-averaging proce- 
dure (from 13 to 19 rotations in our units) used by Fux 
we can compare a particular u-v plot from his figure to 
the time average of all panels in our figure [T] For an 
angle of 3 0° (second row, second column in figure 12 of 
we can see two features stretching from the 
central clump to radial velocities of tt w ±0.2110 ~ ±40 
km/s. Just as expected from our results, the feature 
at positive u (corresponding to negative u in our plots) 
reaches more negative tangential velocities than the one 
moving away from the galactic center. Similar structure 
is apparent for other pattern speeds and angles as well. 
In contrast, in figure 13 of the same paper, where the ve- 
locity distributions are time averaged between 55 and 65 
bar rotations, the features we described above are almost 
not apparent, as the disc is more r elaxed at later times . 

A third example can be found in lAntoja et al.l (|2009| 1 , 
who examined the effect of the Galactic bar and spi- 
ral structure. Similarly to the present study, they did 
not time-average over particle orbits. Although no time 
development of the u-v plane was presented, it is easy 
to identify some structure at low velocities in their Fig. 
1(f). Note also that their discs were unrelaxed, similar 
to those of iMinchev et al.l ()2009( ). which is probably the 
reason for the extended feature (in the u direction) at 
{u,v) « (-20,-30) km/s. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have examined the effect of a central 
bar on the structure in the stellar velocity distribution in 
the solar neighborhood. Unlike in previous works, here 
we do not time average over particle positions but follow 
the time evolution of the system. We find that in ad- 
dition to the central clump in the u-v plane, two strong 
features appear after the formation of the bar, which shift 
with time toward more positive and more negative radial 
velocities (see figure [T]). This is surprising since a system 
subjected to a steady state periodic perturbation must 
not evolve with time due to the presence of an isolating 
integral of motion (the Jacobi's integral). We explain 
these results as the initial response of the disc following 
the bar formation. We associate the streams in the u-v 
plane with two quasi-periodic orbital families precessing 
around the a;i(l) and xi(2) closed orbits near the 2:1 
OLR of the bar. The timescale of this libration is de- 
pendent on the bar's strength as At ~ e~^/^, consistent 
with a simple Hamiltonian model for the resonance (eq. 
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[3]). We find that phase wrapping timescales are longer 
than Hbration timescales thus the disc takes a long time 
to relax near the 2:1 OLR of the bar. 

We have shown for the first time that the Milky Way 
bar could be responsible for low-velocity structure in the 
solar neighborhood velocity field. The stream caused by 
the orbital family just inside the OLR of the bar (fig- 
ure [8]) can be associated with the Pleiades moving group 
seen in the Hipparcos velocity distribution. Similarly, 
the orbital family just outside the OLR can account for 
either the Coma Berenices or the Sirius velocity stream, 
depending on the bar parameters and time since its for- 
mation. We note that the position of the Hercules stream 
(figure [8]) does not change with time as in the case of 
the lower velocity streams. This model requires a bar 
pattern speed of fib/r^o — 1-82 ± 0.07 and a bar orienta- 
tion 30° ^ (/>o ^ 45°. Since the process is recurrent, we 
achieve a good match to the velocity distribution about 
every six rotations of the LSR. However, in order to be 
consistent with the fraction of stars in the Pleiades veloc- 
ity stream, we favor a bar formation time 1.75-2.25 Gyr 
ago. This i s consistent with the observational results by 
iCole andW cinbcrg (2002]), who estimated that the Milky 
Way bar is likely to have formed more r ecently than 3 
Gyr ago. In addition, the recent work bv lMinchev et all 
([2009) in which the authors considered the effect of a 
Milky Way merger on the kinematics of stars in the SN, 
estimated that the Galactic disc was strongly perturbed 
about 2 Gyr ago. Such a merger could have triggered 
bar formation. 

With the upcoming data from the Gaia mission, it will 
become possible to see how the positions of the streams 
we here attribute to the effect of the Milky Way bar, 
change with Galactic azimuth and Galactocentric dis- 
tance, and thus test our model. Furthermore, it might 
be possible to identify the a;i(l) and x\(T) orbital fam- 
ilies in pencil-beam surveys such as ARGOS, BRAVA, 
and APOGEE by observing, for example, in the low 
absorption window at a Galactic longitude of ^ 290° 
(jCarraro and Costal l2009f) . and identifying struct ure in 
the way suggested bv iMinchev and Quillen 1)20081 ). 

This work confirms that resonant features in ve- 
locity space can provide tight constraints on nonax- 
isymmetric disc structure parameters, as shown by 
I Quillen and MinchevI ()2005f ) for the case of spiral density 
waves. We expect that spiral structure can also cause 
precession of resonant orbital families in the reference 
frame of the perturber as seen here, and thus provide 
an insight on spiral structure formation. Preliminary re- 
sults from spiral density wave simulations show that this 



is indeed true (Minchev et al. 2010 , In preparation). 
However, due to the more narrow resonant widths com- 
pared to the case of the bar, this effect might not be as 
important unless the solar radius is almost at the exact 
resonance with spirals. 

Decreasing the pattern speed moves the bar's reso- 
nances outwar ds, causing orbits to be cap t ured into res- 
onan ce (e.g. iRomero-Gomez et al.l l2006t iBaglev et al.l 
|2009() . Therefore, it is likely that the features in the u-v 
plane will remain strong for longer periods of time than 
what is inferred from figure [T] Consequently, the model 
presented here might also be consistent with an old bar, 
provided the bar pattern speed is changing. Further in- 
vestigation exploring this problem is necessary. An addi- 
tional perturbation from the Milky Way spiral structure 
will change the dynamics as well (Quillen and Minche3 
'2005':' De Simone et al."2004:'A ^ia et al.l2009( ). Models 
incorporating both spiral and ba r structure, as the one 
considered bv lChakrabartvl ()2007[) but without suppress- 
ing the effects of the initial disc response, are needed as 
well. We have recently shown that such a configuration 
provides a very powerful stellar radial migration mecha- 
nism (Minchev and Faniaey 2009), however, the effect on 
the u-v plane is yet to be explored. We note that the pa- 
rameter space in the case of two perturbations becomes 
enormous, considering the timescales of each perturbing 
agent associated with the initial response of the disc, in 
addition to the relative orientation and individual time 
evolution of spirals and bar. Thus a thorough investiga- 
tion of the effect of spiral structure (including transient 
spirals) and transient bars is required before a superpo- 
sition of these is considered. Even more complications 
to the disc dynamics will be added if the effect of an 
orbiting satellite is considered (|Quillen et al.ll200"9| ). 

The model in this paper argues against a common 
dynamical origin for the Hyades and Pleiades mov- 
ing groups, in agreement with the results of anal- 
ysis of the Hipparco s catalogue ()Famaev et al.l 120081 : 
iBovv and Hogj l2009l) . This raises the question of 
whether it is the MW spiral structure responsi- 
ble for the Hyades as ex pected from the work by 
I Quillen and MinchevI ()2005f ). A work investigating this 
problem is underway (Minchev et al. 2010, In prepara- 
tion). 
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